Table 1. 

Fraction of Color and Morphology-selected Objects 



I 


Full 


Color 


FWHM 


Combined 


Color/ 


Combined / 


Combined/ 
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Selected 
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Color 
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M31 Halo Field 
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144 


102 


75 


55 


0.71 
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0.38 


18-19 
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136 


89 


64 


0.71 


0.47 


0.33 


19-20 


237 
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104 


66 
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0.51 


0.28 


20-21 


430 


223 


176 


101 


0.52 


0.45 


0.23 


21-22 


805 


379 


314 


184 


0.47 


0.49 


0.23 


22-23 


1384 


629 


592 


333 


0.45 


0.53 


0.24 


23-24 


2156 


992 


1069 


543 


0.46 


0.55 


0.25 
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1996 


1052 


1103 
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0.53 


0.59 


0.31 
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17-18 


64 


47 


37 


30 


0.73 


0.64 


0.47 


18-19 


109 


54 


53 


33 


0.50 


0.61 


0.30 


19-20 


211 


105 


72 


50 


0.50 


0.48 


0.24 


20-21 


298 


127 


103 


65 


0.43 


0.51 


0.22 


21-22 


488 


189 


149 


75 


0.39 


0.40 


0.15 


22-23 


837 


275 


326 


116 


0.33 


0.42 


0.14 


23-24 


1262 


373 


594 


186 


0.30 


0.50 


0.15 


24-25 


880 


287 


442 


135 


0.33 


0.47 


0.15 



Table 2. 

Observed Star Counts and EGM Predictions 



Field 


Area 
[arcmin 2 ] 


X a 

[kpc] 


[kpc] 


Observed 
7=20-22 7=22-24 


EGM Predictions 6 
7=20-22 7=22-24 


This Study 


210.0 


-1.7 


-18.7 


122 c 


543 c 


421 (421) 6141 (6140) 


G302 d 


3.3 


-2.5 


-6.7 


339 


1756 


132 (119) 1955 (1748) 


G312 d 


3.3 


+0.5 


-11.1 


82 


436 


33 (32) 486 (479) 


Gl e 


1.1 


+33.6 


+4.6 


33 


140 


10 (6) 157 (91) 



a (X, Y) = Projected distance of field center from M31's center along its major and minor axes, respectively 
b Numbers in parentheses indicate spheroid stars 

c The star selection efficiency for our 7? = 19 kpc sample is estimated to be about 30% for the brighter apparent 
magnitude range and probably as low as 10% - 20% for fainter stars 

d From the study by Holland et al. (1996) 

e From the study by Rich et al. (1996) 
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ABSTRACT 



Deep UBRI images of a 15' x 15' field in the outer spheroid of the M31 galaxy 
q \ obtained using the Kitt Peak National Observatory 4-meter telescope, and I 

band images obtained with the Keck 10-meter telescope and Low Resolution 
Imaging Spectrograph, are used to isolate a sample of candidate red giant branch 



stars located at a projected radial distance of 19 kpc along the minor axis. 
These stars are distinguished from the more numerous distant field galaxies on 
the basis of broadband U — B, B — R, and R — I colors and image morphology: 
We isolate objects whose colors are consistent with the long, but relatively 
^ ' narrow locus occupied by red giants in UBRI color space (as observed in 

Galactic globular cluster giants and predicted by models spanning a wide range 
of metallicities and ages), and those whose angular sizes are consistent with the 
stellar point spread function, FWHM = 0%-0'!9 (Keck) and FWHM = O'.'9-l '.'5 
(Kitt Peak). We carry out the same analysis of data on a comparison field 
with a similar Galactic latitude to the M31 halo field. The color- magnitude 
diagram of objects in the comparison field is well described by a superposition 
of foreground Galactic dwarf stars (in keeping with a standard empirical model 
of the Galaxy) against a backdrop of contaminating faint blue field galaxies 
(J ^ 21, B — I ~ 1-2.5), while the M31 halo field contains a clear excess of faint 
red objects (I ~ 20 - 23, B — I ~ 2 - 3.5) in addition to these two components. 
The location of this population of faint red objects in the color-magnitude 
diagram is as would be expected for red giant stars at the distance of M31. 
The surface density of red giant candidates in the R = 19 kpc M31 halo field is 
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consistent with the findings of two recent Hubble Space Telescope studies. The 
data indicate that M31's stellar halo is much denser and/or larger than that of 
the Galaxy: (p^f / p^^)(A/1.5)~ v ~ 10, where A is the ratio of the radial scale 
lengths of M31 and the Galaxy and v = —3.8 is the assumed power law index of 
the density profile; in fact, M31's profile may be steeper than this [y < —3.8). 
The color and slope of the red giant branch in M31's outer halo are suggestive of 
a relatively metal-rich population, [Fe/H] ^ —1, in agreement with the Hubble 
Space Telescope measurements. 

Subject headings: galaxies: individual: Andromeda galaxy [Messier 31 (M31), 
NGC 224, UGC 454, CGCG 535-017] - galaxies: formation - stars: red giants - 
stars: metallicity 

1. Introduction 

Studying the metallicity gradient of galactic spheroids is crucial for understanding 
their formation and evolutionary history. The dissipational collapse model (Eggen et al. 
1962; Larson 1974) predicts a strong metallicity gradient because of progressive chemical 
enrichment during the collapse of the protogalactic gas cloud. The accretion model (Searle 
& Zinn 1978), on the other hand, predicts no strong gradient because star formation 
largely precedes assembly of the galactic spheroid. The metallicity gradient of a halo can 
be measured in both the globular cluster system and in field stars; in fact, it is important 
to study both populations independently since they may be dynamically distinct from each 
other. 

Globular clusters are easily observed out to large distances in the Galaxy because of 
their high luminosity and this allows a reliable determination of the metallicity gradient 
in the system of Galactic globular clusters. Zinn (1993) finds evidence for two subsystems 
among clusters: an 'old halo' population which displays a radial [Fe/H] gradient, a small age 
spread, and significant rotation, and a 'young halo' population which displays no metallicity 
gradient, a large age spread, and very little rotation. He associates the former population 
with the initial dissipational collapse of the Galactic halo (a la Eggen et al. 1962) and 
the latter with subsequent accretion of stellar systems (a la Searle & Zinn 1978). On the 
other hand, Carney et al. (1990) find a mean metallicity of [Fe/H] = —1.72 and an absence 
of a strong metallicity gradient as a function of apogalacticon distance over the range 
8 kpc < i? apo < 40 kpc in a proper motion selected survey of 227 Galactic field halo stars 
with velocities lagging the local standard of rest's circular velocity by at least 150 km s -1 . 
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It seems logical to test the universality of the conclusions about galaxy formation 
which are being drawn from the single case of the Galaxy. The Andromeda galaxy (M31) 
provides an external perspective of a large spiral similar to our own and yet is close enough 
for individual stars to be resolved. Huchra et al. (1991) showed that there is evidence for a 
weak metallicity gradient in a sample of 150 M31 globular clusters, with a mean metallicity 
of [Fe/H] = —1.2, which is slightly higher than the mean value of [Fe/H] = —1.4 for Galactic 
globular clusters. 

Over the last decade, several groups have tried to determine the metallicity of field red 
giant branch (RGB) stars in M31's spheroid, starting with the early work of Crotts (1986). 
Mould & Kristian (1986) found a metallicity of [Fe/H] = —0.6 for stars in a field located at 
a projected distance of R = 7 kpc along the minor axis, and assumed that the observed 
color spread was due to a large metallicity dispersion. Mould (1986) failed to determine the 
metallicity at R = 20 kpc, because the contaminating field galaxies could not be removed 
by the standard method of morphological selection and statistical subtraction. Pritchet & 
van den Bergh (1988) estimated [Fe/H] = —1.0 at R = 8.6 kpc on the minor axis, and again 
attributed the large color spread to an intrinsic spread in [Fe/H]. Christian & Heasley's 
(1991) data in an R = 16 kpc field around the globular cluster G219 suggested a large 
metallicity spread, with a mean value [Fe/H] ^ —1.0, as did Davidge's (1993) study at 
R = 6.7 kpc, Durrell et al.'s (1994) study, and Couture et al.'s (1995) study which targeted 
fields around five M31 globular clusters. While these studies have furthered our knowledge 
of M31's halo, sample contamination by distant field galaxies, foreground Galactic dwarfs, 
and M31 disk giants poses a serious obstacle in attempts to determine the metallicity of 
M31's outer halo (R £ 10 kpc). 

Two recent studies of the M31 spheroid have taken advantage of the excellent angular 
resolution of Hubble Space Telescope (HST) Wide Field/Planetary Camera 2 images to 
separate stars from distant field galaxies. In the first study, Rich et al. (1996) analyze 
archival HST images of a field centered on the M31 globular cluster Gl. The line of sight 
to this field (Fig. |l|) intersects M31's highly inclined disk (i ~ 77°; de Vaucouleurs 1958) at 
a relatively large radius -Rdisk ~ R = 40 kpc on the major axis (the value of -Rdisk could be 
different if M31's outer disk is warped). The second study is that of Holland et al. (1996) 
who analyze deep HST images of two fields around M31 globular clusters located 7.6 kpc 
(G302) and 10.8 kpc (G312) from the center roughly along the SE minor axis (Fig. [1]). 
Both studies find the mean metallicity of M31's halo RGB stars to be comparable to that of 
47 Tuc ([Fe/H] = —0.7), with a spread of nearly 2 dex, suggesting that there was a greater 
degree of pre-enrichment during the assembly of M31's spheroid than in the case of the 
Galaxy's spheroid. We will return to a more detailed discussion of these studies later in the 
paper (Sec. 4.2). 



Schematic Plot of M31 Target Fields 
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Fig. 1. — Schematic plot showing the approximate location of the field studied in this paper, 
along with Holland et al.'s (1996) G302 and G312 fields and Rich et al.'s (1996) Gl field. 
The size of each square symbol marking field position corresponds to the actual area of the 
field studied, but not to the field shape. The ellipses indicate the orientation of M31's disk 
(i = 77°, PA = 40°, with finite thickness) and halo (5:3 flattening, same PA as disk); the 
sizes of the ellipses are arbitrary. The 10% change in the cos(<5) projection factor from the 
top to the bottom of the plot has been ignored in this schematic representation. 
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In this paper, we probe RGB stars in the M31 halo using an alternative method for 
discriminating between stars and faint galaxies; one that uses a combination of broadband 
UBRI colors and morphological information (Gould et al. 1992, hereafter referred to as 
GGRF). Early results have been reported by Reitzel et al. (1996a,b). The field of study 
presented here is at a projected distance of R = 19 kpc from M31's center along the minor 
axis (only Rich et al.'s Gl field is at a comparable radial distance out in the spheroid 
assuming it has a 5:3 flattening) and the line of sight to the field intersects M31's disk 
at .Rdisk ~ 75 kpc, further out than in any of the previous studies (Fig. []]). The area of 
the field in this study, 210 arcmin 2 , is more than 60 times larger than each of the Holland 
et al. fields, and more than 40 times larger than the Rich et al. Gl field. Moreover, the 
field stars targeted in this study are well away from any of M31's globular clusters. The 
observations and data reduction procedure are described in Sec. 2. The technique used to 
isolate M31 halo RGB stars from background field galaxies and foreground Galactic dwarfs 
is described in Sec. 3. In Sec. 4, we discuss the properties of M31's spheroid in the context 
of the recent HST studies and outline plans for future work. The main points of the paper 
are summarized in Sec. 5. 



2. The Data 

2.1. Observations 

We present observations of an outer spheroid field in M31 and a comparison field 
carried out with the the Kitt Peak National Observatory (KPNO) 4-meter telescope and a 
2048 x 2048 Tektronix CCD during 4 nights in November 1993. The M31 halo field: 

aj2ooo (M31 halo) = 00 h 48 m 51!0; <W) (M31 halo) = +40°28'05" (1) 

is located 84' from the galaxy's center, at a projected radial distance of 18.8 kpc for an 
assumed distance modulus of (m — M) = 24.43 mag or D = 770 kpc (Freedman & Madore 
1990), at a position angle of 124.8° (E of N) putting it close to M31's SE minor axis (Fig. [[]); 
the position angle of the M31 disk major axis is 40°. The comparison field: 

« J200 o (Comp) = 07 h 49 m 00!0; 5 J200 o (Comp) = +60°03'46" (2) 

provides a control sample: its Galactic coordinates [I 11 (Comp) = 157°, 6 (Comp) = +30°] 
are roughly similar to those of the M31 halo field [7 /7 (M31 halo) = 122°, 
6 /7 (M31 halo) = —22°]. We compare the number, magnitude distribution, and color 
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distribution of Galactic stars in the two fields in the context of the Bahcall & Soneira (1984) 
model (see Sec. 3.4 for details). 

The data set consists of images through U, Bj, R, and I filters with bandpasses 
centered at 3600, 4500, 6500, and 9000 A, respectively. Note, this bandpass system is 
roughly similar, but not identical, to the standard Johnson/Kron/Cousins bandpass system. 
The pixel scale is 0'/47 and the field of view of each individual CCD frame is 16' x 16'. 
The FWHM of the seeing disk ranges from 0'/9 - 1'/5 across the data set, and is generally 
smallest for / band exposures. The primary observations consist of short (15 min in U, 
10 min in each of BjRI), slightly disregistered exposures of the M31 halo field: a total of 
21 exposures in U, 14 in Bj, 13 in R, and 14 in / corresponding to effective total exposure 
times of 5.3 hr in U and 2.2 - 2.3 hr in each of BjRI. The total exposure times are 
somewhat shorter for the comparison field: 15 exposures in U (3.8 hr), 11 in Bj (1.8 hr), 
8 in R (1.3 hr), and 9 in J (1.5 hr). 

In order to empirically verify the locus of RGB stars in multicolor space, we use short 
UBjRI exposures of 6 Galactic globular clusters spanning a range of metallicities. The data 
for the clusters M15, M79, and M2 are of adequate depth and quality; we are unable to use 
the data for the other 3 clusters, NGC 7006, NGC 6934, and Pal 11, as the surface density 
of cluster giants is too low relative to contaminating field Galactic dwarfs. Short UBjRI 
exposures of standard star fields (Landolt 1983) are used for photometric calibration of the 
data. The globular cluster and standard star exposures are from two photometric nights 
during the observing run. In addition, a sequence of zero-second bias frames and short I 
band twilight flat exposures are used for data processing. 

For the purpose of measuring angular sizes of faint objects, we use I band images of 
the M31 halo field obtained using the Keck 10-meter telescope and the Low Resolution 
Imaging Spectrograph (LRIS, Oke et al. 1995) during 1 night in September 1995. The 
4 Keck/LRIS exposures, 1 x 10 min and 3x5 min, form a (mostly) non-overlapping 2x2 
mosaic covering about 75% of the field of view of the KPNO image. The scale of the Keck 
images is 0"22 pixel" 1 and the field of view of each CCD frame is 5'.7 x T.3. The FWHM 
of seeing is % - '.'9 for the Keck observations. Limited use is made of these Keck data; 
unless specifically mentioned, all subsequent data references in the paper pertain to the 
main KPNO data set. 
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2.2. Flat Fielding, Fringe Removal, & Coaddition of the Data 

Each CCD image is overscan subtracted, trimmed, and bias subtracted in the usual 
way. The U, Bj, and R images are then flat fielded using "dark sky flats". Each dark sky 
flat is the median combination, with ±2<r rejection of cosmic rays, stars, and galaxies, of all 
the disregistered images in that band of both M31 halo and comparison fields (20-35 images 
per filter). The / band images are flat fielded differently from those in UBjR as the raw 
CCD frames in / contain Fabry-Perot fringing of night sky emission lines. An / band flat 
field image is created by median combining twilight flat exposures with ±2cr rejection; 
these images are dominated by the (solar) continuum and are largely free of fringing. All / 
exposures are first flat fielded with the twilight flat, which corrects for spatial variation of 
quantum efficiency across the CCD but leaves the fringe pattern intact. A fringe template 
is created by median combining the twilight-flat-fielded I exposures of the M31 halo 
and comparison fields with ±2<r rejection, and by reducing the resulting image to zero 
mean. The amplitude of the fringe pattern is found to vary from image to image (roughly 
0.5%- 1.5% of the sky level), so the fringe template is scaled individually to optimize fringe 
removal for each I band image. 

The images are geometrically transformed to correct for slight image distortion 
(^ 10 -3 ) and aligned for coaddition. Each image is scaled to correct for any variation in 
transparency as a function of airmass and/or time (typically ^ 10%), as determined from 
measurements of secondary photometric standard stars on the image. An additive offset 
is then applied to match the sky background level over common areas of all images in a 
given band, correcting for time variations in the night sky brightness (most extreme in 
the / band). The transparency- and background-matched images in each band are median 
combined with ±2<r rejection of cosmic rays. Cosmic rays are then masked from each 
image by comparing the image to the median image in that bandpass. Finally, all images 
in a given band (UBjRI) are averaged, excluding masked pixels. The combined images 
each cover an area of about 15' x 15'; the dithering between exposures slightly reduces the 
overlap area compared to the field of view of individual CCD frames. 



2.3. Artifact Correction 

The combined CCD images suffer from two kinds of artifacts in the vicinity of bright 
stars (B ^ 17, / ^ 15): charge bleeding and amplifier hysteresis. These problems are most 
severe in the / band image (as the brightest stars are quite red in color and because the 
instrumental throughput is highest in this band) and least in the U band image. While 
the artifacts cover only a small fraction of the total image area (~ 1%), they hamper the 
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performance of the automated object detection and photometry software over a larger 
surrounding area. We apply a set of techniques that reduce, but do not eliminate, the effect 
of artifacts in the M31 halo field image. We do not correct artifacts in the comparison field 
as it is less severely affected than the M31 halo field (higher Galactic latitude, and hence 
fewer bright stars). The affected areas of the comparison field (~ 6% of its area), and 
the artifact-corrected M31 halo field (~ 8% of its area) are properly accounted for in the 
subsequent analysis (Sec. 2.5). 

Extreme saturation of pixels near the center of the image of a bright star can lead to 
the overflow of charge along CCD columns (y axis). We interpolate the sky background row 
by row across the affected region by fitting a linear function in x with a ±2<r rejection to 
avoid stars/galaxies. While this procedure obviously does not recover any of the data in the 
pixels occupied by the bleed, it replaces each bleed with a smoothly varying background, 
thereby allowing faint objects to be detected in the neighborhood of bright stars. 

Amplifier hysteresis results in a temporary change in the effective gain immediately 
after reading out saturated pixels (~ 1% decrease in R and /, <; 10% increase in B and U), 
with the gain gradually returning to its nominal level. This produces prominent streaks of 
low or high values extending for several hundred pixels (2' - 6') along the x axis to the right 
of all bright stars (i.e., downstream in the readout sequence). Multiple streaks caused by 
the numerous bright stars in the M31 halo field combine to form an uneven background and 
affect the probability of faint object detection. Each of these streaks is rectified as follows: 
(1) A linear function in x, f Te {(x), is fit to a 10 pixel-wide (in y) reference region above and 
below the streak, extending along the full length of the streak in x, with ±2<r rejection of 
stars/galaxies; (2) A linear function in x, fi(x), is similarly fit to each row i of the streak; 
and (3) The difference, f re {(x) — fi(x), is added to each row (i) of pixels in the streak. 
This technique does an excellent job of correcting for amplifier hysteresis while preserving: 
(a) the intensity distribution of a star or galaxy whose image happens to lie, partly or 
wholly, within a streak and (b) possible large scale gradients in the sky background level. 

2.4. Detection and Photometry of Stars and Galaxies 

After rectifying charge bleed and amplifier hysteresis artifacts in the data, we construct 
a single deep image for each of the M31 halo and comparison fields. This image is 
constructed from the weighted sum of the images in all 4 bands. The weights are chosen 
to maximize the signal-to-noise ratio for a faint object whose UBjRI colors are roughly 
equal to those of a typical faint RGB star. The software package focas (Jarvis & Tyson 
1981; Valdes 1982) is used to analyze the M31 halo field and comparison field data in the 
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following way: (1) A catalog of individual objects (stars and galaxies) is derived from the 
deep, summed U+Bj+R+I image using a matched filter search technique (with search 
parameters tuned to maximize faint object detection efficiency while minimizing the fraction 
of spurious detections), and a limiting isophote is defined for each object; and (2) For each 
star/galaxy, this limiting isophote is applied to the image in each of the 4 bands (UBjRI) 
to obtain photometry — i.e., color measurements are based on these "common optimized 
apertures". We have checked that, for a set of relatively isolated objects, the focas apparent 
magnitude agrees to within ^ 0.03 mag with the "total" apparent magnitude derived from 
a standard curve of growth analysis. 

After accounting for the difference in line-of-sight reddening between the M31 halo and 
comparison fields, E(B — V) = 0.085 and 0.038, respectively (Burstein & Heiles 1982), 
we find a small offset between the UBRI loci of bright (foreground) stars between the 
two fields corresponding to a residual magnitude zeropoint difference in the B band. This 
difference is probably caused by errors in the transparency correction; empirical checks 
indicate that the difference is not likely to be a result of PSF differences or errors in the 
assumed E(B — V). We apply a slight adjustment to the B magnitude scales of both data 
sets (0.1 mag) to bring them into agreement with each other and with the theoretical 
Bertelli et al. (1994) UBRI stellar loci (Sec. 3.2). This zeropoint adjustment does not affect 
any of the conclusions in this paper — the intrinsic widths of the color selection regions are 
^ 0.1 mag (Fig. [|) and all magnitude measurements are assumed to have an associated 
systematic error of ±0.1 mag in the data analysis (Sec. 2.5). 

We have carried out tests using three different photometry techniques: focas, aperture 
photometry with curve of growth corrections, and PSF-fitting using the digital stellar 
photometry program daophot (Stetson 1987, 1992). The results obtained from these 
techniques are not vastly different. The focas limiting isophotes are better matched to the 
light distribution of galaxies of widely varying shapes and sizes than fixed apertures and 
suffer less from contamination in crowded fields. The focas "common optimized aperture" 
color measurements are somewhat more robust than the daophot measurements for the 
resolved galaxies in the data set; the two techniques yield comparable results for unresolved 
objects. 

The Galactic globular cluster images are analyzed with the help of the daophot 
applying the find and allstar tasks in the usual way. We use a linearly variable point 
spread function (PSF) template to account for slight distortion in image quality across 
the image. Position-dependent aperture corrections are applied to convert the apparent 
magnitudes returned by allstar to total magnitudes. 

The FOCAs-based (M31 halo and comparison fields) and DAOPHOT-based (Galactic 
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globular clusters) total instrumental magnitudes are transformed to the standard 
Johnson/Kron/Cousins UBRI system as defined by the set of Landolt (1983) photometric 
standards. The transformation parameters in each band, derived from observations 
of standard star fields, include a zeropoint accurate to <^ 0.05 mag, a shutter time 
constant of 0.07 s (same for all bands), an independent airmass coefficient for each of the 
two photometric nights, and a linear (Johnson/Kron/Cousins) color term. The color term, 
particularly large for the £?™ str — > ^Johnson convers i ori) j s applied iteratively starting with 
the instrumental color as an initial guess. We hereafter use the term "U BRI magnitudes" 
to refer to calibrated magnitudes on the Johnson/Kron/Cousins system (Landolt 1983). 



2.5. Photometric Accuracy and Completeness 

We estimate the photometric error associated with the measurement of UBRI 
magnitudes of stars/galaxies in the M31 halo and comparison fields. In the M31 halo field, 
the Poisson error is 0.05 mag for objects with U = 23.0, B = 24.1, R = 23.3, and / = 22.6, 
and increases to 0.25 mag for objects with U = 25.5, B = 26.7, R = 25.9, and / = 25.1. 
These estimates are based on shot noise in the photons detected from the object, and on 
the area of the limiting focas isophote and the shot noise per pixel in the sky counts. The 
overall error in the final calibrated magnitude is likely to be somewhat larger than the 
Poisson error because of the following: flat-fielding error, residual image artifacts, error in 
measuring the value of the local sky background, neighbor contamination (varying degrees 
of contamination among the 4 bands can lead to color errors), unaccounted for light outside 
limiting isophote (this does not affect color measurements to first order), and calibration 
error. An empirical check of the measured UBRI colors of Galactic dwarf stars in the 
M31 halo and comparison fields shows that they agree with the Bertelli et al. (1994) model 
predictions to within ~ ±0.05 mag; nevertheless, we conservatively assume that these other 
sources of error contribute an additional (systematic) uncertainty of 0.1 mag (over Poisson 
error) in all magnitude measurements (see Sec. 3.2). 

Based on the sky noise level in the final M31 halo field images, we compute 3<r limiting 
apparent magnitudes within an aperture of area 28.3 pixels or 6.25 arcsec 2 (this is roughly 
equal to the minimum allowable isophotal area in the focas object detection algorithm): 
Un m = 25.6, B\ [m = 26.7, Rn m = 25.7, and I\ im = 25.0. Note, these limiting magnitudes 
are based on a fixed aperture size, whereas the error estimates in the previous paragraph 
are based on the actual isophotal area of each object. The remainder of the data analysis 
presented in this paper is restricted to the sample of stars/galaxies that satisfy the criteria: 
B < B\ im , R < i? lim , and / < I\ im . As discussed in Sec. 3.2 below, the U band is the most 
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critical of the 4 bands in discriminating between faint galaxies and red giants, even in cases 
where there is only an upper limit to the apparent U brightness of an object. We have 
therefore chosen not to impose the U limiting magnitude criterion in defining the sample. 
Instead, any object whose measured U brightness falls formally below the threshold is 
treated as a non-detection in the U band, and is assigned a well-defined upper limit to its 
brightness, U > Un m , for the purpose of color selection (Sec. 3.2). 

Even after applying the above cuts in apparent BRI magnitude, the sample of stars 
and galaxies detected in the M31 halo and comparison fields is not complete over all regions 
of the CCD image: focas's object detection routine has difficulty finding objects in areas 
of the U + B + R + I image surrounding very bright galaxies and stars, especially those 
with charge bleed artifacts. The nominal area covered by the M31 halo data set (field of 
view of KPNO image) is about 227 arcmin 2 but a fraction of this area, (1 — P ) = 0.08, is 
effectively lost. Even though P is probably a weak function of apparent magnitude, it is 
adequate for our purposes to assume it is a constant. We estimate Po by comparing the 
total number of objects detected down to the completeness limits over the entire KPNO 
image to the expected total number based on the observed surface density of objects over 
the same magnitude range in "clean" areas of the image (i.e., areas free of very bright 
stars/galaxies, charge bleeds, and hysteresis artifacts). 

Since the majority of faint objects in the M31 halo field are distant field galaxies 
(an isotropic population), we can estimate the degree of completeness of the sample by 
comparing the observed counts, N(m), to counts of a statistically complete sample of 
galaxies derived from the high-quality, Hubble Deep Field (HDF) data set (Williams et al. 
1996). For the purpose of this comparison only, the above limiting magnitude criteria 
are ignored and the entire KPNO data set is used. The HST instrumental magnitudes 
of HDF galaxies are converted to the Johnson/Kron/Cousins UBRI magnitude system 
adopted in this paper using the transformation relations given by Holtzman et al. (1995); 
the observed distribution of galaxies in the B — R vs B — / plane (in both M31 halo and 
comparison fields) is used to interpolate each HDF galaxy's R band magnitude. For the 
4 bands, the apparent magnitude at which the completeness drops to half the maximum 
value is: U50 = 26.3, -B50 = 27.6, R50 = 26.5, and /50 = 25.8. These values are significantly 
fainter than the 3a limiting apparent magnitudes at which the data set is truncated; the 
completeness fraction is essentially 100% for objects brighter than the limiting magnitudes. 
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2.6. Angular Size Measurements 

The angular diameter (#fwhm) of each object detected in the M31 halo field and 
comparison field is measured using the final, coadded KPNO / band images, as they have 
the best seeing of the 4 bands (FWHM of seeing is about 1 ".2 and 1 "5 in the two fields, 
respectively). The seeing in a few of the individual 10 min KPNO / exposures is as good as 
0'/8 (FWHM), but these images do not have adequate signal-to- noise for the majority of the 
objects which are quite faint. For each object, 6*fwhm is determined by fitting a Gaussian 
to the azimuthally-averaged radial intensity profile. The outer parts of the stellar PSF are 
not very well described by a Gaussian so the #fwhm values reported in this paper may be 
systematically different from measurements based on a more realistic PSF template; this 
study however relies only on the relative angular sizes of objects (Sec. 3.3). 

We also determine widths of the best fit Gaussian profiles for objects that lie in the 
area of the M31 halo field covered by the Keck/LRIS / band images. The Keck data set 
has better angular resolution (0 '.'6 - '.'9 seeing FWHM and '.'22 pixel" 1 scale) than the 
coadded KPNO image and is of comparable or greater depth (the factor of 6 larger primary 
mirror area and better seeing more than compensate for the much shorter exposure times). 
We use Keck-based #fwhm measurements in preference to those derived from the KPNO 
data whenever available; the former allow improved morphological distinction between stars 
and compact galaxies at faint apparent magnitudes. 

3. Searching for M31 Halo Red Giants: Needles in a Haystack 

3.1. Defining the Stellar Locus in UBRI space 

Accurate definition of the loci of RGB stars in U BRI space is very important for 
the color selection technique used in this paper. Following the method used by GGRF 
who determined stellar loci using the Yale stellar evolution models (Green et al. 1987), we 
determine RGB stellar loci using the Bertelli et al. (1994) stellar evolution models. The 
Bertelli et al. model isochrones match the color distribution of Landolt (1983) standard 
stars extremely well, in the Johnson/ Kron/ Cousins UBVRI photometric system defined 
by these stars. Since the photometric calibration of the KPNO data set used in this 
paper is tied to these same Landolt standards (in terms of magnitude zeropoints and color 
transformations), the Bertelli et al. models are a natural choice for defining the stellar loci 
used in this study. 

We use Bertelli et al. model isochrones bracketing the 5-20 Gyr range of ages and 
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Fig. 2. — The loci of stars in the BRI and U BR color-color diagrams. [All colors are adjusted 
to the amount of reddening appropriate for the M31 halo field, E(B — V) — 0.085.] Left 
panels — Theoretical isochrones from Bertelli et al. (1994) for ages of 5-20 Gyr and for 
metallicities of [Fe/H] = —1.70 to [Fe/H] = (bold dots). The dashed line in the lower left 
panel corresponds to a 1 Gyr stellar population with super-solar metallicity, [Fe/H] = +0.4. 
The isochrones exclude main sequence stars with M ^ 0.6 M causing an artificial truncation 
(e.g., at B — R — 2.4, isochrones U — B = 1.2). The solid lines define the "selection 
parallelogram", the region used to select RGB stars in the M31 halo and comparison fields. 
The reddening vector is shown for A v = 1 mag (its location in the color-color plane is 
arbitrary). Right panels — Distribution of relatively metal poor stars in the Galactic globular 
clusters, M2 (squares), M79 (triangles), and M15 (circles), along with the RGB selection 
regions shown in the left panels. The distribution of cluster RGB stars is in good agreement 
with theoretical models; the feature extending beyond the blue end of the selection region 
(B — R < 0.8) is formed by cluster horizontal branch stars. 
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for the metallicity range [Fe/H] = —1.70 to [Fe/H] = (bold dots in left panels of Fig. ^) 
to define the stellar selection loci in the R — I vs B — R and U — B vs B — R color-color 
planes. The isochrones are reddened by an amount appropriate for the M31 halo field, 
E(B — V) — 0.085 (Burstein & Heiles 1982), using a standard Galactic interstellar dust 
extinction law (Cardelli et al. 1989). The selection region (the narrow region in Fig. | with 
parallel sides which we refer to as a "parallelogram" for simplicity) is defined as the region 
that completely encompasses all of the RGB points in the isochrones. While we do not 
expect there to be many young or high metallicity objects, we have chosen the selection 
region to avoid biasing the resulting sample to any particular metallicity or age. 

Even though the stellar locus is very extended in any one color (e.g., it spans the range 
0.8 < B — R < 3.3), the different colors are strongly correlated and this makes the locus 
narrow in color-color space. Changes in metallicity and age tend to move stars mostly 
along the length of the stellar loci: this is especially true in the BRI color-color diagram; 
increasing [Fe/H] does make the U — B color slightly redder at fixed B — R, so the width of 
the UBR color-color selection region is correspondingly wider to account for this. 

We perform two independent empirical checks to verify the accuracy of the stellar 
locus definitions and the photometric conversion of instrumental magnitudes to the 
Johnson/Kron/Cousins system: 

(1) Globular cluster data, obtained in the same UBjRI system in conjunction with 
the M31 halo field and comparison field observations, yield a distribution of RGB 
stars that is in good agreement with the Bertelli et al. (1994) model isochrones. 
This is shown in the right panels of Fig. ^] (the sequence of stars beyond the left 
edge of the selection box corresponds to the cluster horizontal branch). Prior to the 
comparison, the colors of the globular cluster stars are dereddened and then adjusted 
to the line-of-sight reddening estimated for the M31 halo field, EiB — V) = 0.085 
(Burstein & Heiles 1982). Note however that data are only available for the relatively 
metal-poor clusters, M15, M2, and M79, thereby restricting the check to the range: 
-2.5 £ [Fe/H] £ -1.5. 

(2) As discussed in Sec. 3.2 below, the distribution of foreground Galactic metal-rich 
([Fe/H] ~ 0) disk dwarf stars in UBRI color-color diagrams of the M31 halo and 
comparison fields is consistent with Bertelli et al. isochrones of solar or super-solar 
metallicity (compare dashed line in lower left panel of Fig. ^| to lower left panel of 
Fig. This check is indirect however because the foreground stars are (mostly) main 
sequence stars whereas the selection loci are optimized for red giants. 
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3.2. Color Selection 

The main source of contamination for the study of M31's faint spheroid stars is the 
population of background faint field galaxies. These galaxies are expected to be about as 
common as M31 halo stars at faint magnitudes (there are 25 arcmin -2 down to B = 26) 
and a fair fraction have small angular sizes, which makes it difficult to distinguish them 
from stars on the basis of image morphology alone. The expected surface density of M31 
spheroid red giants in this 19 kpc field is less than 26 arcmin -2 down to / = 24 (B ~ 26) 
as discussed in Sec. 4.2. GGRF have developed a method which uses accurate UBRI CCD 
photometry for isolating faint stars (B < 26) and rejects most of the background galaxies 
on the basis of colors while eliminating only a fraction of the stars. The essence of the 
method lies in the fact that stars occupy a relatively tight locus in any color-color diagram 
(Fig. |2|), while galaxies exhibit a broad distribution of colors. 

Figure |3] demonstrates the basis for the color selection method via BRI and UBR 
color-color plots of all objects detected in the M31 halo field. For an object to appear in 
these color-color plots and in the AU vs AI color excess plots below (Fig. f|) it must satisfy 
the U < Un m criterion along with the usual BRI criteria; this additional restriction causes 
there to be substantially fewer objects in these plots than in the full sample (cf. Figs. |6] 
and |7|), especially at the faint end of the distribution. The left panels include all objects 
brighter than the expected tip of the M31 RGB, / < 20.5. Most of the bright objects in this 
low latitude = 22°) M31 halo field are foreground Galactic stars as evidenced from the 
fact that they lie close to the defined selection region. The selection regions are, however, 
optimized for red giant stars of solar or lower metallicity, while most of these bright objects 
are disk dwarfs with high (possibly super-solar) metallicity. Thus it is not surprising that 
the I < 20.5 objects deviate from the selection region slightly, especially in the UBR 
diagram. The right panels of Fig. ^| show all objects with / > 20.5, which includes any M31 
RGB stars as well as a large number of distant field galaxies. It is obvious that the galaxies 
occupy a more extensive region of the plot than do the stars. In the BRI diagram, the 
galaxies are slightly redder in R — I (<~ 0.4 mag) than the stellar locus at a given B — R 
color; the separation is clearer in the UBR diagram: for example, the typical galaxy is 
^ 1 mag bluer in. U — B than the stellar locus for B — R ^> 2. 

Following GGRF we use each object's measured B — R color and the stellar locus to 
predict the colors, (R — I) pTe d and (U — -B) pre d, an RGB star would have if it had the same 
B — R color as the object. The measured colors of the object are then used to calculate 
color excesses: 



AI = -[(R-I)-(R-I) pred \ 



(3) 
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AU=(U-B)-(U- B) pred (4) 

The color excesses measure the vertical distance in the BRI and UBR color-color diagrams 
of each object from center line of the stellar locus "parallelogram" (Fig. |2]). The U data are 
pivotal in the application of the color selection technique to red giants, even though the 
low intrinsic stellar brightness and relative inefficiency of the instrument in the ultraviolet 
conspire to make U band detection of distant RGB stars a challenging prospect. As 
described in Sec. 2.5, objects that are fainter than the limiting magnitude in U (but 
which satisfy the apparent brightness criteria in BRI) are assigned a U band color excess 
AU > AU lim , where Af/ lim = (C/ lim - B) - (U - B) pred . 

Figure |4] shows the distribution of color excesses. The relatively bright objects 
(17 < I < 20.5) in the M31 halo field shown in the left panel lie near but not exactly at 
the origin. This is only to be expected: while most of these bright objects are foreground 
Galactic stars, they are main sequence stars and not red giants. In the center and right 
panels, we plot the fainter objects (20.5 < I < 25) in the M31 halo and comparison fields, 
respectively. Note the broad and asymmetric distribution of these faint objects in (AI, AU) 
space most of which belong to the faint background field galaxy population (especially in 
the comparison field, right panel). 

Color selection is accomplished by requiring objects to have colors consistent with 
both BRI and UBR stellar loci. An object is considered to pass color selection only if 
its color excesses (AJ, AU) are such that the 1.5 a measurement error ellipse intersects a 
rectangular reference region centered on the origin (Fig. |4j), where the width and height 
of the rectangle are equal to the vertical width of the "parallelogram" selection regions in 
the BRI and UBR color-color plots, respectively (Fig. |2|). In other words, an object is 
assigned a color-based classification of "star-like" only if 5 < 1.5, where S is the significance, 
in units of the photometric error, of the separation of the object from the point closest to 
it within the central rectangular area in (A/, AU) space. Like GGRF, we determine the 
measurement errors for AI and AU jointly, as there is a strong covariance between these 
quantities. For each U band non-detection that is brighter than the limiting magnitudes 
in BRI, we check whether the combination of its AJ value and AU > A U^ m criterion is 
consistent with the rectangle at the origin. 

3.3. Morphological Selection 

The UBRI color selection described above yields a sample of objects consisting of M31 
RGB stars, foreground Galactic stars (mostly disk dwarfs), and any field galaxies whose 
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Fig. 3. — Color-color diagrams for the M31 halo field: R — I vs B — R (upper panels) 
and U — B vs B — R (lower panels). Point size is scaled with apparent I magnitude, with 
larger dots indicating brighter objects. The left panels show objects with 17 < I < 20.5. 
These bright objects are mostly foreground main sequence stars in the Galactic disk; their 
distribution matches the selection region (see dashed line in Fig. |2[), though not perfectly as 
selection is optimized for lower metallicity RGB stars. The right panels show objects with 
20.5 < I < 25. The majority of the objects in this faint magnitude range are galaxies which 
have a much larger spread in color than the stars. Note, most galaxies with B — R ^ 2 are 
~ 1 mag bluer in U — B than the selection region. 



-18- 




Reduced color: M Reduced color: AT Reduced color: M 



Fig. 4.— The distribution of color excess indices, AI and AU. Dot size scales with 
apparent / magnitude: large dots = bright objects, and vice versa. The left panel shows 
objects with 17 < I < 20.5 in the M31 halo field: most of these are Galactic dwarfs with 
high metallicities, with colors that deviate slightly from the expected locus of RGB colors. 
The center and right panels show objects with 20.5 < / < 25 in the M31 halo field and 
the comparison field, respectively; the majority of these are faint galaxies (especially in the 
comparison field) exhibiting a broad distribution of color excesses. The width and height 
of the bold rectangle at the origin correspond to the vertical widths of the selection regions 
used in the BRI and U BR planes, respectively. 
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colors happen to be consistent with stellar loci. We also apply the more traditional method 
of star-galaxy separation based on image morphology. 

Figure [| shows the distribution of angular sizes (^fwhm) of objects in the M31 halo 
field as measured on the KPNO image and on one of the Keck/LRIS I band images. About 
half the objects in the sample are well resolved relative to the stellar PSF, #fwhm > #crit 
(morphologically "galaxy-like"); the rest are assigned a morphological classification of 
"star-like". We adopt 9 cvit = 1 '.'36 for the KPNO image, and 6 crit = 0(73, 0'.'84, 0'.'86, and 
0'/88 for the 4 Keck/LRIS images. These 8 cr i t values are slightly larger than the FWHM of 
the PSF in the center of the corresponding CCD image to allow for variations in the PSF 
across the field of view: the stellar image quality is degraded near the corners of the field 
due to optical distortion/focus variation. For an object larger than a few arcseconds, the 
measured size is typically inaccurate since it is based on a Gaussian fit. This does not affect 
our analysis, however, since 6*fwhm for such objects is invariably large; hence the objects 
are correctly classified as morphologically galaxy-like. 

Keck/LRIS-based star-galaxy separation is used in preference to KPNO-based 
star-galaxy separation in areas where both measurements are available (75% of the M31 
halo field); the superior angular resolution and signal-to- noise ratio of the Keck J-band 
data result in improved star-galaxy discrimination, especially in the case of compact, 
barely resolved galaxies and blended groups of faint objects. It is possible to compare 
star-galaxy discrimination based on Keck vs KPNO #fwhm measurements in the area where 
the two data sets overlap: the morphological classification is star-like by both methods for 
34% of the objects, galaxy-like by both methods for 26%, star-like according to the Keck 
measurement and galaxy-like according to the KPNO measurement for 27%, and vice versa 
for 13%. 

The distribution of measured angular sizes may be used to investigate whether the 
color selection technique is working. In the unselected sample, 56.0% of all objects with 
/ < 20 have #fwhm values that lead to their being classified as morphologically star-like, 
whereas after color selection 61.6% of the objects are morphologically star-like. For 
/ < 23, only 44.4% of all objects have stellar profiles, yet 53.4% of the color selected 
objects are star-like. For / < 25, 48.8% of the full (non-color selected) sample has star-like 
morphology, compared to 55.3% for the color selected sample. The fact that the fraction 
of morphologically star-like objects is higher in the color selected sample than in the full 
sample is indicative that color selection is indeed picking out stars. The effect is subtle 
since neither morphological selection nor color selection alone is perfectly efficient. A 
fraction of the stars (~ 30%) happen to be scattered beyond 5 > 1.5 due to photometric 
error, and there are compact galaxies that fail the UBRI selection criterion as well as 
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Fig. 5. — The distribution of angular sizes for objects detected in the M31 halo field as 
measured on the KPNO image (upper panel) and on one of the four Keck images (lower 
panel). The threshold for morphological selection, 6> crit , is indicated by a vertical line. The 
solid curve shows the entire (unselected) sample; the dotted curve is for objects whose colors 
are consistent with those of an RGB star (5 < 1.5), while the dashed curve is for the rest of 
the objects (5 > 1.5). Note, objects with 5 < 1.5 tend to be more compact than those with 
5 > 1.5, an indication that the color selection does preferentially select stars. 
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extended galaxies whose colors happen to intersect the RGB locus. Furthermore, errors in 
the measurement of angular size (due to object/sky photon noise, neighbor contamination, 
etc.) cause #fwhm for some stars to scatter beyond crit , especially at faint magnitudes. 

There is nothing special about our use of #fwhm measurements for star-galaxy 
discrimination nor about our particular choices of cr i t values. After carrying out tests of 
focas's Bayesian star-galaxy separation scheme on the data sets, we choose to adopt the 
simpler #fwhm measurement scheme; the effect of noise on star-galaxy discrimination is 
easier to quantify in the latter scheme. It is worth noting however that results based on the 
focas morphological classification are not qualitatively different from those based on #fwhm 
measurements. The 6*FWHM-selected sample no doubt includes some compact galaxies for 
which the measured angular size is less than the critical value. Conversely, a small fraction 
of stars are scattered across to #fwhm > #crit due to measurement error in the angular size 
caused by noise and/or the presence of neighbors. We have chosen crit values to achieve a 
balance between contamination vs incompleteness. 



3.4. Detection of Red Giants in M31 

The detection of M31 RGB stars against the numerous background field galaxies and 
foreground Galactic stars is difficult, but the combined power of color and FWHM selection 
yields a sample of M31 halo RGB stars. The pre-selection color-magnitude diagram (CMD) 
of all objects detected in the M31 halo field is shown in Fig. |6|(a); it is practically impossible 
to see any concentration of M31 RGB stars in this plot. After FWHM selection a hint of 
concentration of M31 RGB stars becomes visible in the CMD [Fig. 0(b)]. Morphological 
selection removes 51.2% of the objects in the full sample, most of them relatively well 
resolved galaxies, yet contamination from compact galaxies and foreground main sequence 
stars remains. The UBRI selection technique removes 49.7% of the full set of objects 
[Fig. H(c)]. Most foreground lower main sequence stars with B — I > 3 are removed (these 
deviate systematically from the RGB loci) as are a significant fraction of the distant field 
galaxies (these have a broad distribution in color-color space), yet galaxies whose UBRI 
colors happen to be close to those of a RGB star survive selection. While a concentration 
of M31 RGB stars in the CMD is again revealed, as for the FWHM selected sample, the 
exact location of the concentration is uncertain due to contamination. When FWHM and 
color selection are applied simultaneously, the complimentary nature of these selection 
screens reveals a fairly clean sample of M31 RGB stars [Fig. §](d)]. Comparing the Galactic 
globular cluster RGB fiducials of M5 and 47 Tuc [plotted in Fig. |6|(b) and (c), but not 
in (d) so as not to guide the reader's eye] to the main concentration in the FWHM- and 
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?7_Bi?/-selected CMD [Fig. |(d)], it is clear that the objects detected are in a location of 
the CMD appropriate for M31 RGB stars. 

We carry out a detailed comparison between the M31 halo field sample and the sample 
of objects found in the comparison field. The final image of the comparison field has an area 
of 254 arcmin 2 ; its effective area is 94.1% of this, or 239 arcmin 2 (P has been estimated in 
the same way as for the M31 halo field — see Sec. 2.5). The M31 halo field has a smaller 
effective area (88%) than the comparison field. The latter data set is scaled down in order 
to match the effective area of the M31 halo field by removing a random 12% of the objects 
from the sample. This reduces the number of objects in the full (unselected) comparison 
field sample to 4217. The matched comparison field data set is then processed through the 
same UBRI color selection and FWHM selection screens as the M31 halo field data set. 
The comparison field data were obtained under slightly worse seeing conditions than the 
M31 halo field data; the KPNO / band comparison field data set is used for morphological 
selection with 9 cr \ t = l'/65. No Keck/LRIS images of the comparison field are available for 
FWHM measurements. 

Figure ^ shows I vs B — I CMDs of the M31 halo and comparison fields (upper and 
lower panels, respectively), clearly demonstrating the presence of an additional component 
of red giants in the former data set. The left panels [(a) and (b)] contain the full data 
sets prior to any selection. Even in this non-selected sample, it is obvious that the M31 
halo field contains many more objects than the comparison field, 7469 versus 4217. This is 
mostly due to the additional population of M31 RGB stars, although as we discuss below, 
the M31 halo field does contain a larger number of foreground Galactic stars since it is 
at a somewhat lower Galactic latitude. In the center panels [(c) and (d)] we compare the 
fields after UBRI selection: 50.3% of the M31 halo field objects survive, while only 36.1% 
of the comparison field objects survive, an indication of the larger population of objects 
with star-like colors in the former field. The additional population of M31 RGB stars is 
plainly visible, even though many galaxies with star-like colors remain. The right panels 
[(e) and (f)] contain only objects that have both star-like colors and morphology. There is 
a noticeable additional component of M31 halo stars in Fig. |7|(e). As much as 57.0% of the 
FWHM selected M31 halo field objects survive the additional color selection, while only 
41.1% survive in the comparison field. 

The objects that survive color and morphological selection in the comparison field 
[Fig. |7j(f)] represent a mix of Galactic stars and compact galaxies that happen to have 
star-like colors. RGB fiducials for the Galactic globular clusters M5 and 47 Tuc are also 
plotted to indicate the region of the CMD where M31 RGB stars are expected to be 
concentrated. The significant concentration of additional objects in the M31 halo field 
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Fig. 6. — (a) Color-magnitude diagram of all M31 halo field objects. (b) CMD after 
morphological selection, (c) CMD after U BRI color selection, (d) CMD after U BRI and 
FWHM selection. The number of objects in each panel is indicated in the upper left. Most 
of the bright objects (I ^ 20) are foreground Galactic stars, the prominent vertical edge at 
B — I = 1 representing the main sequence turnoff, while the majority of faint objects are 
distant galaxies. The diagonal concentration of faint red objects that survive in (d) resembles 
the RGB fiducials of Galactic globular clusters displaced to Dm3i = 770 kpc shown in (b) 
and (c): [L->R] M5 ([Fe/H] = -1.4) and 47 Tuc ([Fe/H] = -0.70). 
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Fig. 7. — (a+b) Color-magnitude diagrams of all objects in the M31 halo and comparison 
fields, respectively, with the latter scaled down to correct for the difference in the effective 
areas of the two data sets. (c+d) Same as (a+b) after UBRI selection. (e+f) Same 
as (a+b) after UBRI and morphological selection. The number of objects in each panel 
is indicated in the upper left. The unselected M31 halo field contains many more objects 
than the corresponding comparison field sample but it is difficult to identify the additional 
population in the former sample [compare (a) vs (b)]. Note the significant concentration 
in (e) of additional objects (RGB stars) in the final UBRI- and FWHM-selected M31 halo 
field data set. This may be compared to the RGB fiducials of M5 and 47 Tuc plotted in (f), 
as for Fig. ||. 
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sample [Fig. 0(e)] is indeed at this expected CMD location; the M31 halo field also contains 
the foreground and background populations that are present in the comparison field CMD. 
The UBRI- and FWHM-selected M31 halo field sample contains 1321 more objects than 
the corresponding comparison field sample. The number of these additional M31 halo 
field objects increases towards fainter apparent I magnitudes (see Table 1). There are 
145 additional objects in the range 20 < / < 22, and 574 additional objects in the range 
22 < I < 24. Even in the magnitude range brighter than the expected tip of the M31 RGB, 
/ < 20, there are 116 more objects in the M31 halo field than in the comparison field. This 
is due to the higher surface density of foreground Galactic stars in the M31 halo field. 

The degree of foreground contamination due to Galactic stars in the M31 halo 
and comparison fields is investigated using an empirical model of the Galaxy (IASG 
model — Ratnatunga & Bahcall 1985). The model incorporates standard density profile 
parameters and metallicity distributions for the Galaxy's spheroid, thin disk, and thick 
disk to compute the predicted stellar color and apparent magnitude distributions. These 
model CMDs are compared to the CMDs of the M31 halo and comparison fields. While the 
IASG model has not been tested at faint magnitudes, the bright end stellar distributions 
agree well with the M31 halo and comparison field data [Figs. ||(a), (b), (c), and (d)]. Note 
the "wall" of main sequence turnoff stars near the left edge of each panel. The predicted 
density of stars with / < 20 and B — I < 3.5 agrees with that observed in the corresponding 
fields. The IASG model seems to over predict the number of stars with B — I > 3.5; in fact, 
recent star count analyses by Gould et al. (1996, 1997) indicate that the abundance of low 
luminosity red (B — I > 3.5) main sequence stars in the Galaxy is lower by a factor of ~ 3 
than that predicted by the IASG model. Another reason for the discrepancy between the 
M31 halo CMD and the IASG model CMD [Figs. |(a) vs (c)] is that many of these faint 
red stars are screened out of the M31 halo sample during UBRI selection because the locus 
of such stars departs from the RGB locus on which the selection "parallelogram" is based 
(see left panels of Figs. ^ and |3|). 



3.5. Statistical Removal of Contaminants from M31 RGB Sample 

We use the color- and morphology-selected, area-matched comparison field data set 
[Fig. |8](b)] to assess the amount of contamination in the corresponding M31 halo field 
sample. Each object in the comparison field CMD is matched to the object "closest" to it in 
color-magnitude location in the scaled M31 halo CMD, provided it is located within a local 
elliptical region. This elliptical region is arbitrarily defined by semi axes lengths of 0.3 mag 
in B — I color and 0.5 mag in apparent I magnitude. The matched pairs are removed 
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from both CMDs. Only 81 objects remain unmatched out of the 757 in the UBRI- and 
FWHM-selected comparison field sample (~ 10%), while 1402 out of 2078 (67%) remain 
unmatched in the scaled M31 halo field. This preliminary sample of "excess" objects in the 
M31 halo field shows a strong concentration around the expected CMD location of RGB 
stars at the distance of M31 [Fig. |(e)]: CMD|e = CMD^ - CMD|b. 

The M31 halo field is closer to the Galactic plane (|6| = 22°) than the comparison 
field (|6| = 30°), and is expected to contain 1.7 times more foreground Galactic dwarf stars 
than the area matched comparison field sample on the basis of the IASG model [compare 
Figs. H(c) and (d)]. We carry out another statistical subtraction between the two IASG 
model CMDs to estimate the number and distribution in color-magnitude space of the 
extra foreground Galactic stars in the scaled M31 halo field. We also simulate the fall off 
in UBRI selection efficiency for low luminosity, red main sequence stars by excluding an 
increasing fraction of stars with B — I > 3.5; the excluded stars are indicated by crosses in 
Figs. |8](c) and (d). This sample of extra IASG model stars is then used to (statistically) 
remove residual foreground contaminants from the preliminary statistically-subtracted 
M31 halo field sample to obtain the final sample of M31 halo RGB candidates [Fig. |8](f)]: 
CMDgf = CMD|e - (CMD| - CMD|d). 

3.6. Star Selection Efficiency 

The efficiency of the UBRI- and FWHM-based star selection procedure may be 
estimated by studying the fraction of objects selected at every step of the procedure, 
coupled with a priori knowledge of the ratio of stars to galaxies, as a function of apparent 
brightness. Table 1 shows the selected numbers and fractions of objects for both M31 
halo and comparison fields in 1 mag bins in the / band. At bright magnitudes (/ < 19), 
the surface density of field galaxies is only 0.43 arcmin 2 (Lilly et al. 1995a,b) or about 
90 galaxies in the KPNO field, and foreground Galactic stars are expected to outnumber 
field galaxies 7:2 in the low latitude M31 halo field. Color selection appears to retain about 
two-thirds of all stars in the brightest magnitude bins. This is the expected fraction enclosed 
within the 1.5er ellipse (see UBRI selection criterion in Sec. 3.2) for a two dimensional 
Gaussian distribution. Morphological selection picks out ^ 50% of the remaining objects. 
This fraction is significantly less than unity (even at bright magnitudes where the effect of 
galaxy contamination is negligible) because of measurement error in the angular size (due 
to photon noise, crowding, etc.) coupled with the fact that the PSF quality is not constant 
over the CCD image. While increasing the value of cr i t would include more stars, it would 
also correspondingly increase the number of barely resolved galaxies that pass FWHM 
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Fig. 8. — Color-magnitude diagrams illustrating the statistical removal of contaminants. 
The number of objects in each panel is indicated in the upper left. RGB fiducials are the 
same as in Fig. ||. (a+b) Same as for Fig. |7|(e+f). (c+d) CMDs derived from IASG model 
predictions for the M31 halo and comparison fields, respectively. The faint crosses indicate 
the fraction of faint red lower main sequence stars that fail the color selection test; their 
number is given in parentheses, (e) Statistical excess of objects in the M31 halo field CMD 
in (a) over that in (b). (f) Same as (e) corrected for excess foreground contamination in 
the M31 halo field [(c) minus (d)]. 
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Table 1. 

Fraction of Color and Morphology-selected Objects 



I 


Full 


Color 


FWHM ' 


Combined 


Color/ 


Combined/ 


Combined/ 




Sample 


Selected 


Selected 


Selection 


Full 


Color 


Full 


M31 Halo Field 


17-18 


144 


102 


75 


55 


0.71 


0.54 


0.38 


18-19 


192 


136 


89 


64 


0.71 


0.47 


0.33 


19-20 


237 


130 


104 


66 


0.55 


0.51 


0.28 


20-21 


430 


223 


176 


101 


0.52 


0.45 


0.23 


21-22 


805 


379 


314 


184 


0.47 


0.49 


0.23 


22-23 


1384 


629 


592 


333 


0.45 


0.53 


0.24 


23-24 


2156 


992 


1069 


543 


0.46 


0.55 


0.25 


24-25 


1996 


1052 


1103 


620 


0.53 


0.59 


0.31 


(Matched) Comparison Field 


17-18 


64 


47 


37 


30 


0.73 


0.64 


0.47 


18-19 


109 


54 


53 


33 


0.50 


0.61 


0.30 


19-20 


211 


105 


72 


50 


0.50 


0.48 


0.24 


20-21 


298 


127 


103 


65 


0.43 


0.51 


0.22 


21-22 


488 


189 


149 


75 


0.39 


0.40 


0.15 


22-23 


837 


275 


326 


116 


0.33 


0.42 


0.14 


23-24 


1262 


373 


594 


186 


0.30 


0.50 


0.15 


24-25 


880 


287 


442 


135 


0.33 


0.47 


0.15 
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selection. Overall the combined selection efficiency for stars at bright magnitudes is ~ 30%. 

At magnitudes fainter than about / = 20, the UBRI selection step picks only ~ 50% of 
all objects. It should be noted though that this is not a direct measure of the star selection 
efficiency. The measured fraction of selected objects is decreased (relative to the 'true' 
UBRI selection efficiency for stars) by the presence of field galaxies in the unselected sample 
which dominate the sample at the faintest magnitudes. Some of these galaxies happen to 
have star-like colors, however, and these contaminants tend to increase the color-selected 
fraction. The increased importance of systematic, non-Gaussian errors at faint magnitudes 
(e.g. due to interference from close neighbors, sky background measurement error) likely 
causes more than a third of the stars to escape color selection. FWHM selection at faint 
magnitudes retains ~ 50% of the color-selected subsample, but this is not necessarily a good 
measure of the morphological selection efficiency for stars for reasons explained above. The 
true FWHM selection efficiency at faint magnitudes is expected to be lower than 50% (the 
estimated efficiency at bright magnitudes) because of increased error in the measurement of 
#fwhm- Thus, the overall selection efficiency for faint stars is likely in the range 10% -20%. 

The comparison field data are also useful for investigating the selection efficiency. As 
noted earlier, this field is at a slightly higher Galactic latitude than the M31 halo field 
and therefore contains a smaller number of foreground Galactic dwarf stars. Thus the 
star-to-galaxy ratio is lower in the comparison field even for magnitudes brighter than 
/ = 20, a range in which there is no contribution from M31 red giant stars in the M31 halo 
field sample. It is therefore no surprise that the fraction of color- and morphology-selected 
objects is somewhat lower in the comparison field than in the M31 halo field over the full 
range of apparent brightnesses (Table 1). A more detailed comparison between the selection 
efficiencies in the M31 halo and comparison fields is precluded by subtle differences in their 
faint end noise characteristics and in the seeing FWHM (which affect the UBRI- and 
FWHM-selected fractions, respectively). 

4. Discussion 

4.1. Estimating the Metallicity of M31's Outer Spheroid 

The metallicity of the population of M31 halo stars detected is estimated by comparing 
the color-magnitude distribution of the U BRI- and FWHM-selected, statistically-subtracted 
sample of objects to the RGB fiducials for the Galactic globular clusters M5 and 47 Tuc 
in Fig. ||(f). Clearly the fiducial for M5 is not a good match as most of the M31 halo 
field objects brighter than / = 23 are redder, and hence more metal rich, than M5 
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([Fe/H] = -1.4). The clump of faint (I > 23), blue {B - I ~ 1 - 2) objects in this plot 
is produced at least in part by residual contamination by distant faint blue field galaxies 
[compare Figs. |(c) and (f)]. The fiducial for 47 Tuc is a better match to the data, and 
yet there is a substantial spread about this fiducial. A substantial fraction of the objects 
with I < 23 lie to the right of the 47 Tuc fiducial which indicates that the mean metallicity 
of M31's outer halo is comparable to, and possibly a bit higher than, that of 47 Tuc 
([Fe/H] = —0.70), with a fairly large spread in color. 

The mean photometric errors estimated from Poisson noise are ob-i — 0.02 for objects 
near the tip of the M31 RGB (/ ~ 20.5) and <Jb-i — 0.18 near the limit of the sample 
(/ ~ 24). The observed spread in color is too large to be accounted for simply in terms 
of photometric error. Despite our efforts to remove all possible contaminants, some of the 
observed spread could be due to residual contamination from foreground Galactic stars and 
background faint compact galaxies. The IASG model has not been thoroughly tested at 
faint magnitudes; moreover, Poisson fluctuations in the field galaxy counts cause differences 
between the comparison and M31 halo fields. These possible contaminants cannot be 
entirely responsible for the observed color spread however as there is no reason why the 
contaminants should contribute primarily near the observed CMD location of the M31 
RGB. In fact, compact galaxies contribute primarily at faint magnitudes (/ ^ 22) and 
blue colors (B — I < 2) as can be seen by comparing the UBRI- and FWHM-selected 
comparison field sample in Fig. |S|(b) with the IASG model prediction for this field in 
Fig. |8|(d). The uncertainty of the IASG model at faint magnitudes is harder to quantify; 
some Galactic stars may still be contaminating even the foreground subtracted CMD in 
Fig. ||(f). We conclude that much of the spread observed in the RGB concentration is due 
to a true spread in the metallicity of M31 halo stars. Future followup spectroscopy will test 
this conclusion (Sec. 4.3). 



4.2. Comparison of this Study to Recent HST Studies 

In this section, we compare the final, statistically-subtracted sample of RGB candidates 
in the M31 halo CMD to those found in the HST studies of Holland et al. (1996) and 
Rich et al. (1996). Fig. |1] shows the locations of all the fields in relation to M31's disk and 
halo (schematic representation). We adopt the authors' plausible hypothesis that the vast 
majority of objects seen in Rich et al.'s R = 40 kpc (Gl) CMD and in Holland et al.'s 
R = 7.6 kpc (G302) and 10.8 kpc (G312) CMDs are M31 field halo stars (see details below). 
We use the External Galaxy Model (EGM) developed by Hodder (1995) to estimate the 
number of halo and disk stars expected to be seen in each of these fields. The EGM is an 
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adaptation of the Bahcall & Soneira (1984) Galaxy model designed to give star counts and 
color distributions for any field in an external spiral galaxy using a de Vaucouleurs profile 
for the halo and an exponential disk. 

There is a large number of free parameters in the EGM (disk and halo scale lengths, 
density ratio, halo axis ratio, and overall normalization), but a plausible model based on 
47 Tuc's red giant luminosity function and color distribution, a disk exponential scale 
length of 5.5 kpc, a halo effective radius of 2.6 kpc, and an axis ratio for the halo of 0.6 
can be normalized to match the star counts observed in the HST fields to within ~3%. 
Furthermore, the fraction of M31 disk stars predicted by the EGM for the two Holland 
et al. fields agrees with the upper limits to the fractional disk contamination estimated 
by the authors: 0.1 for G302 and 0.03 for G312. The halo effective radius of 2.6 kpc for 
the major axis corresponds to a minor axis effective radius of 1.6 kpc for the flattened 
halo. This is consistent with the minor axis effective radius of 1.3 kpc found by Pritchet & 
van den Bergh (1994). The EGM predicts about 6500 RGB stars with 20 < I < 24 within 
the 210 arcmin 2 effective area of the KPNO image of the M31 halo. Most of these are 
expected to be halo giants with negligible contamination from M31's disk (0.01%). Table 2 
presents the observed number of stars and EGM predictions for this study and for the three 
HST fields. 

The surface density of M31 stars predicted by the EGM in the R = 19 kpc field we 
have studied is lower than in any of the three HST fields, even Rich et al.'s Gl field at 
R = 40 kpc. A 5:3 flattened halo contributes about equally to the Gl field and ours, but 
M31's disk contributes an additional 70% (40% of total) in the Gl field. Contrary to the 
impression one may get from Fig. [I], the degree of disk contamination is much higher in Gl 
than in the inner Holland et al. (G302) field. Even though the G302 field is expected to 
contain more disk stars than Gl, this is more than compensated by the much larger number 
of halo stars in G302, making the fractional disk contamination in G302 quite small. This 
is simply a reflection of the fact that galaxy halos have steeper density profiles than the 
typical disk exponential law. The M31 disk RGB fraction in Gl could be even higher if the 
halo is less flattened (b/a > 0.6) and/or if a thick disk is present. Thus, the large spread 
in the M31 halo metallicity inferred by Rich et al. may in fact be due primarily to disk 
contamination (possibly including stars with super-solar [Fe/H]). 

After FWHM and UBRI selection and statistical subtraction, we detect 665 objects 
down to a limit of / < 24 in the M31 halo field [see Fig. |5](f)], whereas the EGM predicts 
about 6500. This would seem to indicate a ~ 10% selection efficiency, but the efficiency is 
higher than this over most of the magnitude range. About ~ 30% of the number predicted 
by the EGM is observed in the range 20 < I < 22, which is close to estimated efficiency 
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Table 2. 

Observed Star Counts and EGM Predictions 



Field 


Area 


X a 


Y a 


Observed 


EGM Predictions 6 




[arcmin 2 ] 


[kpc] 


[kpc] 


7=20-22 


7=22-24 


7=20-22 7=22-24 


This Study 


210.0 


-1.7 


-18.7 


122 c 


543 c 


421 (421) 6141 (6140) 


G302 d 


3.3 


-2.5 


-6.7 


339 


1756 


132 (119) 1955 (1748) 


G312 d 


3.3 


+0.5 


-11.1 


82 


436 


33 (32) 486 (479) 


Gl e 


1.1 


+33.6 


+4.6 


33 


140 


10 (6) 157 (91) 



a (X, Y) = Projected distance of field center from M31's center along its major and 
minor axes, respectively 

b Numbers in parentheses indicate spheroid stars 



c The star selection efficiency for our R = 19 kpc sample is estimated to be about 30% 
for the brighter apparent magnitude range and probably as low as 10%-20% for fainter 
stars 

d From the study by Holland et al. (1996) 
e From the study by Rich et al. (1996) 



-33- 



of the star selection process at relatively bright magnitudes (as discussed in Sec. 3.6). At 
fainter magnitudes the selection efficiency drops for a variety of reasons: non-gaussian 
photometric errors (effect of crowding on isophotal definition, sky measurement errors, etc.) 
and increased FWHM measurement error (causing stars to scatter beyond CT it)- Overall, 
the number of objects seen in our field is roughly consistent with the EGM model prediction 
given the selection efficiency. 

The estimated metallicity of M31's outer halo from our study (Sec. 4.1) is consistent 
with that found in the recent HST studies. Holland et al. (1996) find a spread in metallicity 
of —2 <^ [Fe/H] <^ —0.2 with the majority of stars having [Fe/H]~ —0.6. The contamination 
due to M31 disk stars is estimated using the EGM to be ~ 0.1 for the G302 field and ~ 0.03 
for the G312 field. This, along with the estimated foreground Galactic contamination of 
only 15 disk stars, cannot completely explain the population of very red (high metallicity) 
objects observed in the G302 and G312 fields. If the disk of M31 is warped, however, the 
fractional disk contamination could be higher in these Holland et al. fields. Rich et al. 
(1996) find the field population around Gl to be as metal rich as 47 Tuc, but warn that the 
apparent high metallicity could also be due to the sample containing objects of intermediate 
age. As discussed above, contamination by M31's disk is expected to be significant in the 
Gl field. 



4.3. Density and Size of the M31 Halo 

The observed RGB counts in the Holland et al. (1996) HST minor axis fields and 
in our R = 19 kpc minor axis field can be used to constrain the physical parameters of 
M31's halo. The halo is assumed to be oblate and viewed edge on — i.e., with its symmetry 
axis in the plane of the sky. A (deprojected) power law index of v = —3.8 fits the surface 
densities observed in the R = 7.6 kpc (G302) and 10.8 kpc (G312) Holland et al. fields 
derived from stars with 20 < / < 23.5. This slope is a bit shallower than but consistent 
with the range of v = — 4 to —5 found by Pritchet & van den Bergh (1994) over the (minor 
axis) radial range of R = 3 - 20 kpc based on stars within 2 mag of the tip of the RGB. 
Pritchet & van den Bergh find that M31's halo density profile steepens with increasing 
radius (in a log-log plot), so it is not surprising that a single power law fit to the inner halo 
(R < 10.8 kpc) is shallower than one that includes the outer halo. 

For a power law density profile, the column density at a projected radial distance R 
along the minor axis is given by: 



34 



r 

N = j_ 



oc 



x 2 + (R/c) 

Rq 



dx 



(5) 



N = n 




(6) 



where c = 0.6 is the halo axial ratio, x is the distance along the line of sight, and 



G(u) 



T(-u/2 - l/2)r(l/2) 
T(-u/2) 



(7) 



The volume density no, at a radial distance Ro in the equatorial plane of the halo, may thus 
be inferred from the observed surface density: 



The halo volume density is calculated in the apparent magnitude range 20 < / < 23.5 
to avoid red clump stars and to directly compare to Morrison's (1993) measurement of the 
local Galactic halo stellar density. A projected density of N — 7870 kpc~ 2 is derived from 
counts of stars (1290) in this / magnitude range in the inner Holland et al. field (G302). 
The area of this field (two Wide Field Camera CCDs) is 3.3 arcmin 2 which corresponds to 
0.17 kpc 2 at the distance of M31. The inner (7.6 kpc) Holland et al. field, G302, is not 
exactly on the minor axis; the ellipse with c = 0.6 that intersects this field has a semi-minor 
axis length of R = 6.8 kpc, and this is the value used in the above equations. 

As M31 is about A 1.5 times larger than the Galaxy, we choose R = 12 kpc in 
order to calculate the halo density at a position in M31 comparable to the solar location in 
the Galaxy. The derived M31 halo density at Rq = 12 kpc is tiq = 341 kpc -3 . Thus, uq 
is much larger than the local Galactic halo density of 36 kpc~ 3 for stars with My < 0.5 
(Morrison 1993). This indicates that M31's halo is much denser, and/or larger than that 
of the Galaxy, (p^IF / Pmw)(^/^-^) u ~ 9.5 ± 0.3, where A is the ratio of the characteristic 
radial scale lengths of M31 and the Galaxy, v ~ —3.8 is the power law density profile slope, 
and the (formal) error bar is based on Poisson error in the star counts in the inner Holland 
et al. field, G302. The uncertainty in the overdensity estimate is substantially higher than 
the formal error bar quoted above when one takes into account the possible range of profile 
slopes, v (see below). 

The actual density of M31's inner halo is likely to be even higher than this estimate: 
We have adopted a single power law profile in the above integral, with a slope of v = —3.8 




(8) 
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based on M31's inner halo, but the outer profile is known to steepen both from the counts 
in our R = 19 kpc field and from Pritchet & van den Bergh (1994); this leads to an 
overestimate of the column density for a given spatial density n , or to an underestimate of 
no since we normalize to the observed column density At any rate, it is clear than the halo 
of M31 is significantly denser/larger than the halo of the Galaxy. 

The density profile of the outer halo of M31 (beyond a minor axis radius of 10 kpc) has 
a slope of v « — 5, as indicated by the EGM de Vaucouleurs spheroid fit to the observed 
counts at R = 10.8 kpc (G312, Holland et al. 1996) and R = 19 kpc (this study), and by 
Pritchet & van den Bergh's (1994) study. This is steeper than the profile of outer Galactic 
halo which may be approximated by a power law index in the range %f = —3.0 to —3.5 
(Zinn 1985; Preston et al. 1991; Gould et al. 1998). There are significant uncertainties 
though in these measurements of the profile of the Galaxy's outer stellar halo, and the 
halo of M31 is not well fit by a power law but rather by a de Vaucouleurs law (Pritchet & 
van den Bergh 1994), so directly comparing the power law indices may be misleading. 



4.4. Future Work/Followup Spectroscopy 

Although it is becoming increasingly clear that the halo of M31 is fairly metal rich, 
an accurate mean metallicity and spread remain to be determined. Contamination from 
foreground Galactic stars and background field galaxies, as well as from disk RGB stars 
in M31, is still non- negligible. The combined color- and morphological-selection technique 
described in this paper can be used effectively to study M31's inner halo where the fractional 
contamination is expected to be lower than in our R = 19 kpc field. High angular resolution 
imaging (e.g., with HST) can help reduce background galaxy contamination. 

We are currently undertaking followup Keck/LRIS spectroscopy of the brightest 
decade of M31 halo RGB candidates in the UBRI- and FWHM-selected sample [Fig. §(d)] 
(Guhathakurta & Reitzel 1998). Candidates for spectroscopic observation are selected 
without regard to color, so as not to introduce any metallicity bias. We do not select objects 
from the statistically-subtracted version of the CMD as it is not complete. Spectroscopy 
should allow kinematic confirmation of the identity of M31 RGB stars: galaxies are 
eliminated altogether on the basis of their redshift and M31 stars are quite well separated 
from foreground Galactic disk stars (t>M3i ~ —300 km s -1 vs vuw disk ~ km s -1 ). A broad 
versus narrow distribution of radial velocities will distinguish M31 halo stars from disk 
stars on M31's minor axis, respectively. In addition to removal of residual contaminants, 
spectroscopy allows a direct measurement of the metallicity of each star using the near 
infrared Ca n line strengths, Wc a - A well developed empirical calibration method exists 
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to determine [Fe/H] from W Ca and the brightness of the star SI relative to the horizontal 
branch level (Olszewski et al. 1991; Armandroff & Da Costa 1991; Suntzeff & Kraft 1996). 

5. Summary 

We analyze deep U BRI images obtained with the KPNO 4-meter telescope of a field 
in the outer halo of M31 (R = 19 kpc in projection on the minor axis of M31), and of 
a comparison field located at similar Galactic coordinates. Keck/LRIS / band images 
of the M31 halo field are used for the purpose of making angular size measurements 
and morphological star-galaxy separation. We use a UBRI color selection technique in 
conjunction with morphological selection to isolate a sample of M31 RGB stars from the 
numerous background population of faint field galaxies. A detailed comparison shows 
that the M31 halo field contains an excess of objects over the populations present in the 
comparison field. An empirical Galactic (IASG) model is used to estimate the foreground 
contamination from Galactic stars in both fields; the predicted number of bright stars agrees 
well with the observed distributions in the M31 halo field and the comparison field. We 
use the comparison field data and the IASG model to statistically subtract contaminants 
from the M31 halo field, revealing a distribution of objects at the expected location of the 
RGB at the distance of M31. These stars have a mean metallicity of [Fe/H] > —0.7 with 
a significant spread about this value, ranging from M5 ([Fe/H] = —1.4) to more metal rich 
than 47 Tuc ([Fe/H] = —0.70). This result is in good agreement with recent HST studies of 
the halo of M31 by Holland et al. (1996) and Rich et al. (1996). The spatial density of M31's 
halo red giants is an order of magnitude greater than the red giant density at a comparable 
radius in the Milky Way halo. Alternatively, the halo of M31 has a characteristic size scale 
almost 3 times larger than that of the Galactic halo. The outer stellar halo of M31 has a 
steeper slope {v ~ —5) than that of the Galaxy [y = —3 to —3.5). 
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